/// 


A)--  n 


X77-1J/37/501 


J 


To  be  PreTOnted  at  the 

Eleventh  Annual  ^ilomar  Conference 
On  Circuits,  Systems,  and  Computers 

Pacific  Grove,  California 

7-9  November  1977 


THE  DISTRIBUTION  OF  THIS  REPORT 
IS  UNLIMITED 


Etectronict  Reiaarch 

Rockwell  International 

3370  Miralomt  Avanu* 

Anahaim.  California  92S03 


V 

D DC 

|R|r?f75f?m  nf?| 

ur^ IAN  90  WW 


EEiEinns 


037/501  EQUIVALENT  GAIN  FORMS  FOR  CTD  FILTERS 


EQUIVALENT  GAIN  FORMS  FOR  CTO  FILTERS 

Dr.  Stanley  A.  White 

Electronics  Research  Center 
Rockwell  International 
Anaheim,  California 

Abstract 

^here  are  three  equivalent  forms  for  expressing  the  gains  of  transversal 
analog  filters  which  are  mechanized  with  chajrge-transfer  devices  (CTO's) 
Due  to  charge-transfer  inefficiency  (CTD.f'n,  the  equivalence  of 
expression  is  not  obvious.  This  paper  presents  a simple  tabulation 
equivalences,  their  uses,  and  their  derivations.  This  problem  and  sol 
also  extend  to  include  the  recursive  case.  Arbitrary  discrete-time 
recursive  transfer  functions  are  often  mechanized  digitally  as  second 
sections.  Equivalent  analog  realization  using  charge-transfer  devices  can 
not  use  the  same  gain  parameters  becuase  of  the  effects  (again)  of  charge- 
transfer  inefficiency.  This  paper  further  includes  an  expression  for 
corrected  gains  to  compensate  for  this  charge-transfer  inefficiency  in  the 
recursive  structure.  ^ 
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A wealth  of  design  material  exists  to  guide 
the  designer  of  discrete-time  transversal  filters.* 
These  design  aids  determine  the  tap  weights,  or 
gains,  for  filters  of  the  form. 

^ -k 

g{z)  ' r a.z  (1) 

k=0 

where  the  gain  values  are  given  by  the  a. 's  and 

-1  ^ 
where  z is  the  unit  delay  operator.  Digital 

implementation  of  such  filters  is  a straight- 

2 

forward  procedure,  but  analog  implementation  is 

3 

more  difficult.  A digital  delay  of  a single 

sample  time,  z*^  is  realized  precisely.  An  analog 
delay  which  uses  charge-transfer  devices  (CTD)  , 
does  not  exhibit  the  ideal  transfer  function  of  z’  , 
but  instead  the  transfer  function  is 

D (z)  . — 

' 1-nz 

4 

where  n is  the  charge-transfer  inefficiency.  If 
one  were  to  attempt  to  mechanize  the  transfer 
function  (1)  by  CTD's  directly  using  the  unmodi- 
fied tap  weights,  he  would  obtain  a new  transfer 
function: 

f(z)  - z a z'*'  (2) 

k*0  ''^l-nz  / 

which  provides  the  desired  transfer  function  only 
when  n'O.  Now,  notice  that  aj^  is  the  gain  which 

is  set  by  actual  multipliers,  but  the  effect  of  the 
n on  the  transfer  function  can  be  considerable  and 
need  be  accounted  for.- 


while  (1)  has  none.  We  have  two  options:  elimin- 
ate the  denominator  in  (2)  by  using  a power-series 
expansion  and  truncating  after  the  Kill  term,  or 
acknowledge  that  the  denominator  is  there.  We 
choose  the  latter  procedure  which  leads  to: 


f(z) 


h.Z  . 

j J 


(z-n)*^ 


9(z) 

(z-n)*^ 


(3a) 


or,  equivalently. 


f(z) 


t c.t'^ 

TTTk 


(1-nz  ')''  (1-nz  ') 


(3b) 


We  are  going  to  design  a CTD  filter  by  first  deter- 
mining the  b's  (or  the  c's)  then  the  actual  multi- 
plier gains  (the  a's).  These  b's  and  c's  are 
obviously  not  the  actual  multiplier  gains;  however, 
these  expressions  of  (3)  and  (4)  are  much  closer  to 
the  "standard"  form  of  (1)  because  their  numerators 
are  exactly  in  the  form  of  (1). 

The  design  procedure  is  now  executed  in  3 
steps: 

1.  The  desired  transfer  function,  f(z),  is 
determined.  We  then  modify  it  to  form  g in  the 
step  below. 

2.  By  establishing  FIR  filter  design  procedures, 
such  as  described  in  references  1 and  2,  the 
gains  b^  or  Cj  are  obtained  from 

g(z)  » f(z)(z-n)'^  » I b.z-^  • g(z) 
j-0  J 


We  wish  to  express  (2)  in  a form  similar  to 
(1)  because  of  the  ready  accessibility  of  design 
aids,  however,  (2)  has  a denominator  of  order  K 
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(4a)  . 


a a 


or 


II.  Derivations 


g(z)  • f(z)(l-nz'^)*^  =■  l c.z"'^  (4b) 

j=0  J 

3.  The  actual  multiplier  gains  a. , are  deter- 
mined from  the  bj  or  c^: 


K 

.. . iu 


K+k 


(1-n)' 


(5a) 


-K+k 


(1-n)" 

which  is  derived  in  Section  II. 


(5b) 


Above  we  have  a synthesis  procedure.  Now 
let's  examine  the  solution  to  the  inverse  problem, 
the  analysis  procedure;  i.e.,  given  the  values  of 
aj^,  determine  the  filter  impulse  and  frequency 

responses.  Existing  software  abounds  to  develop 
the  answer  from  the  b's  or  c's. 

The  analysis  procedure  is  executed  in  4 steps: 

1.  The  filter  structure  is  given  in  the  form  of 

(1). 

2.  The  coefficients  b.  or  c.  are  obtained  from 

the  by  J J 

bj  . j/K-."-'"''' C)  '‘‘I 


and 


"K-j 


b,.  I 
•’  k.K-j  ^ " 
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* b 


K-j 


(6b) 


3.  FIR  filter  analysis  programs  can  be  used  to 
evaluate  the  performance  of: 


g{z)  » r b.z 
•0  J 


j 


J-0 


or 


-K  1 

z g(z)  - I c.z^ 

j-0  J 


(7a) 


(7b) 


4.  Performance  of  the  acutal  filter  can  be  evalu- 
ated from 

f{z)  • g(z)(z-n)*''-  z‘''g(z)(l-nz*’)‘'^  (7c) 


Since  each  individual  linearized  CTD  cell  can 
be  represented  by  the  transfer  function 


n,(z)  = Z-’  - ^ 

’ 1-nz^ 


(8) 


,.th 


the  N order  tapped  delay  line  can  be  described 

by  / \2 


* V*1  ^"2 


'H 


(9a) 

(9b) 


z a,.(1-n)'‘(z-n)'^''' 

■ K K 

(z-n)*^  (z-n)"^ 


The  numerator  of  (9c)  can  be  rewritten  as 


g(z)  = z a.  (l-n)'‘(z-n)'^’'‘ 
k=0  " 


= z a^.^(1-n)'^''‘(z-n)'' 


(10) 


k=0 


by  redefining  k as  K-k.  Next  we'll  expand  (z-n)", 
so  that  now 


g(z)  - z 
k*0 


K-k 


K . K 
r z^  E ak  B.. 
j-0  k=j  J" 


where  j*k-r,  and 


*k  * *K-k^’*''^ 


K-k 


(11) 


(12a) 


and 


“jk 


Equations  (9c),  (11),  (12a)  and  (12b)  combine  to 
yield  (3a)  with  b.  defined  as  in  (6a).  When  (3a) 
K ^ 

is  divided  bv  z , (3b)  and  (6b)  emerge.  Equations 
(10)  and  (7a)  can  be  combined  to  give 

K K 

g(z)  - E a„  ,(1-n)'^*'*(z-n)'^  - E b-z-* 
j-0  j-0  J 

(13) 


then 


1 im  3 


(13) 


Jl K J-k 


yielding 


Tt^  '’j'’ 

J*k 


b.r 


Vk  = j,  (t-n)'^-'' 

J •*  K 


(14) 


F(2^ 


1+Az**^  + Bz'^ 


this  can  be  used  as  a second-order  filter  multi- 
plexed across  N data  streams,  but  cursed  with 
crosstalk.  This  can  also  be  used  as  a single 
channel  filter  with  a periodic  frequency  response. 
In  either  case,  the  following  argument  holds. 


Suppose  that  the  transfer  function  of  an 
ideal  discrete-time  processor  is 


G(z)  . EjiJ- 

' ' ^z) 


(15) 


its  poles  and  zeros  are  given  by  the  respective 
solutions  to 

q(z)  • 0 and  p(z)  > 0 . 

The  singularities  in  each  case  are  given  by 
expressions  of  the  form 


z • z 


k • 


(16) 


If  the  processor  is  to  be  time-shared  between  N 
channels,  we  could  rewrite  (6)  as 


(17) 


Using  the  CTD  mechanization 
iN 


r T n 

v:*  = 

1 -n  . 


(18) 


and  we  see  that  the  filter  singularities  have  been 
modified  to 

iN 


(19) 


We  would  like  to  prewarp  Z|^  to  z!  so  that  (17)  is 
preserved  in  (19),  i.e.,  so  that^ 


^l'  ' ^k 


r i/N  If 

= [(l-n)z;  .nj 


(20) 


Solving  the  above  for  z'  yields  the  prewarped 
singularity  location  ^ 


which  can  be  expressed  in  the  form  of  (5a). 

From  (3a)  and  (3b)  one  can  see  that  b.=  c„  . 

J ^ J 

which  when  substituted  in  (5a)  give  (5b). 

III.  The  Recursive  Case 

Here  we  shall  introduce  one  additional  level 
of  complexity.  For  reasons  of  economics  we  often 
find  CTO  transversal  filters  of  large  order;  CTO's 
come  that  way.  Since  stability  problems  mitigate 
against  large-order  recursive  filters,  the  order  is 
restricted  to  a small  number,  like  two.  In  order 
to  capitalize  on  the  large  number  of  CTO's  which 
one  can  economically  make,  we  find  structures  like 


fz  „ 1 

^k 

. 1-n 


N 


(2i; 


Now  let's  relate  this  to  filter  gain  values.  The 
conventional  description  of  singularities  is 
usually  given  in  polar  coordinates  so  that 


(22) 


Conjugate-complex  pairs  of  singularities,  z and 
zj  , appear  in  the  characteristic  equation:k 


0 = (z-z,^)(z-zjj)  » z ♦ A,^z  + 


where 


A|^  = 2r|^cos0|^  = -2  Re[Z|^] 


and 


(23) 


(24) 


(25) 


and  where  Aj^  and  are  filter  gain  values. 

The  results  of  the  prewarping  will  be  incorporated 
into  the  gain  values: 


2Re 


[^k]  ‘ '^yfk  <=« 

/ -1  ^k  N ^ \ 

IN  tan IC-  1 


and 


r/'" 


7n 


cos 


nr 


(i-n)2 


(27) 
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which  indeed  places  the  singularities  of  the  CTD- 
mechanized  filter  in  their  proper  locations.  For 
large  N we  may  simplify  (26)  and  (27)  by  the 
approximations: 
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9C 


COS 


l-pr  -1/N 


•Vk 


(28) 


and 


^ \ 


1 nr 

l-nr^- 


2N 


(29) 


The  gain  for  real  singularity  is  given  directly 
by  (21). 
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